Exploring multi-stability in semiconductor ring lasers: theory and experiment 
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We report the first experimental observation of multi-stable states in a single-longitudinal mode 
semiconductor ring laser. We show how the operation of the device can be steered to either mono- 
stable, bistable or multi-stable dynamical regimes in a controlled way. We observe that the dynam- 
ical regimes are organized in well reproducible sequences that match the bifurcation diagrams of a 
two-dimensional model. By analyzing the phase space in this model, we predict how the stochastic 
transitions between multi-stable states take place and confirm it experimentally. 



PACS numbers: 42.65.Pc, 42.55.Px,42.60.Mi 

Multi-stability is a general feature of nonlinear systems 
which attracts attention in a broad set of subjects includ- 
ing hydrodynamics [JJ, plasma physics [2], biology [3], 
neural networks [UH], chemical reactions [6] and optical 
systems [7HH]- The phase space of a multi-stable system 
is in general very intricate due to the strongly interwoven 
basins of attractions of the coexisting stable structures, 
and is often further complicated by the presence of struc- 
tures such as chaotic saddles [TOj HTJ- For this reason, 
the dynamics of a multi-stable system is characterized 
by a larger complexity than their bistable counterpart, 
leading to phenomena such as attractor hopping [TTJ [T2] 
or chaotic itineracy [5] . While being of broad interdisci- 
plinary interest, multi-stability is especially interesting in 
the case of semiconductor lasers, due to their large num- 
ber of applications and their wealth of dynamical regimes 
(see [TJ] and reference therein). However, the fast time 
scales involved, the presence of spontaneous emission of 
photons blurring off the structures, the difficulty to con- 
trol the internal parameters, the inaccessibility of some 
dynamical variables, all make the experimental recon- 
struction of the phase space of semiconductor lasers an 
extremely challenging task. Therefore, despite the large 
number of theoretical work [TTJ [THlfTflj , the dynamical 
complexity of multi-stable semiconductor lasers remains 
experimentally unaddressed. 

fn this contribution, we experimentally address the phase 
space of Semiconductor Ring Lasers (SRLs), which are a 
novel class of semiconductor lasers characterized by cir- 
cular geometry [17j. We have focused our investigation 
on SRLs for several reasons. From the theoretical point 
of view, many dynamical regimes including multi-stable 
ones have been predicted [TB] but not observed yet . From 
the technological point of view, an understanding of the 
phase space of SRLs would be highly desirable to improve 
performances in applications such as all-optical memories 
[T§] and allow for a better control of the device [T51 120) . 
We start our investigation with a comprehensive model 
2 1 J , which allows us to predict bifurcation sequences and 




FIG. 1. Experimental SRL set-up. Four waveguide contacts 
are depicted of which only the one denoted by I w is biased. 

the different possible phase portraits of SRL-dynamics, 
including multi-stable ones. Experimentally, we show 
how to control the parameters of a real SRL fully ex- 
ploring its parameter space and the corresponding phase 
space portraits. 

Consider a SRL operating in a single longitudinal 
mode. Two directional modes with intensities P lj2 and 
phases (f>i,2 can propagate in the ring cavity. For sake of 
clarity, we assume that Pi is the clockwise (CW) prop- 
agating mode and Pi the counterclockwise one (CCW) 
[see Fig. flj. The coupling between CW and CCW is 
described by a complex coupling parameter with ampli- 
tude K and phase (j>k [T7]- On time scales slower than 
the relaxation oscillations, the total power P1 + P2 is con- 
served and the dynamics of the SRL can be described by 
the following two-dimensional asymptotic model [21 j : 

= — 2 sin 4>k sin ip + 2 cos <f>k cos yj sin 

+ J sin 9 cos 9, (1) 
cos Oyh = aJ sin 9 cos 9 + 2cos(f>k sin ip 

+ 2 sin <pk cos ijj sin 9. (2) 

where 9 = 2 arctan yJP 2 /Pi - tt/2 G [-tt/2, tt/2] rep- 
resents the partitioning of power between modes, and 
tp — cf>2 — 4>i & [0, 27r] is the phase difference between 
the counter-propagating modes. J is the rescaled bias 
current and a is the lincwidth-cnhanccment factor. A 
full bifurcation analysis of Eqs. ([lJ-Q has previously re- 
vealed a wide range of dynamical states [IS]. We focus 
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FIG. 2. Bifurcation diagrams of Eqs. |T])-([2| for (a) (j> k = 0.2 
and (b) 4>k = 1-5 Stable solutions of Eqs. (Q-|2| are marked 
in black, unstable ones in gray. Dashed lines are used to indi- 
cate periodic solutions. The Roman numbers I - V indicate 
different dynamical regimes. 
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FIG. 3. Phase space portraits of Eqs. |TJ-([2]| for different 
values of J and (j>k- The notation is as defined in the text, 
(a) <t> k = 0.2, J = 1.05, (b) fc = 1.5, J = 0.75, (c) cj> k = 1.5, 
J = 0.5, (d) <f> k = 1.5, J = 0.66. 



here on two different values of the parameter fa represen- 
tative for yet unexplored sequences of dynamical regimes 
including multi-stable ones. 

Consider <j) k = 0.2. The bifurcation diagram for such 
value of fa is shown in Fig. [2ja) while representative ex- 
amples of its phase space are shown in Fig(3] For small 
values of the bias current, the system operates in bidi- 
rectional regime and the phase space of the system (not 
shown) consists of a stable state B\ coexisting with an 
unstable one B2 (Region I) . An increase of the bias cur- 
rent J leads to the appearance of two more stable states 
corresponding to unidirectional CW and CCW rotating 
solutions U\ 2- The SRL therefore operates in a tristable 
regime (Region II). The corresponding phase space is 
shown in Fig. [3ja). Three stationary states U\, U2 and 
B\ coexist. An unstable state B2 corresponding to bidi- 
rectional in-phase lasing and two saddles Si ,2 are also 
present in the system. The basins of attraction of the 
three states are separated by the stable manifolds of Si 
and S2 in such a way that the basin of attraction of B\ lies 
in between the basins of attraction of U\ and t/2- When 
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FIG. 4. Simulations vs. experiments for a sequential transi- 
tion Ui — > Bi — *■ U2 for the phase space described in Fig.|3^b). 
(a) projection on the phase space; (b) numerical time series; 
(c) experimental data. The model parameters for (a)-(b) are 
J — 1.05, cj>k = 1.02, a = 3.5. The experimental conditions 
for (c) were I p = 45.31mA and I w = 9.81mA. 



spontaneous emission noise is introduced in the system, 
spontaneous attractor hopping may appear. In the limit 
of vanishing noise intensity, the topology of the phase 
space predicts only transitions between a unidirectional 
mode and the bidirectional mode. Direct transitions be- 
tween Ui and U2 are possible only for larger values of the 
noise intensity. An example of a simulated hopping tra- 
jectory is shown in Fig. Qa)-(b), and will be compared 
with experimental switches later. Increasing the value of 
the bias current J, the saddles Si. 2 migrate towards B\ 
making the basin of attraction of Bi shrink. Therefore 
the residence time in the bidirectional lasing mode is ex- 
pected to decrease and the laser to operate most of its 
time in cither one of the unidirectional modes. When the 
bias current is increased above a critical value, the sad- 
dles collide with B\ and the bidirectional out-of-phase 
mode become unstable. The SRL is then operating in a 
bistable regime with two stable unidirectional modes U\.2 
(Region III). The corresponding phase space is shown in 
Fig. |3jb). Here U\ and U2 are the unidirectional modes 
whereas B\ is a saddle point. The stable manifold of 
Bi separates the basins of attraction of TJ\ and L^- In 
this regimes, noise-induced hopping is expected between 
the two unidirectional modes, whereas we do not expect 
to observe any residence in the bidirectional regime. A 
further increase of the current leads to longer residence 
times in U12 but no further bifurcations are expected. 

When <j) k = 1.5 the bifurcation curve as shown in 
Fig.[2jh) is qualitatively different from the previous case. 
For small values of the bias current, the system oper- 
ates in the bidirectional regime, similar to Region I for 
4>k — 0.2. Increasing the current above a critical value, 
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the bidirectional operation loses its stability and the SRL 
exhibits periodic oscillations known as alternate oscilla- 
tions 17] between CW and CCW modes (Region IV). 
In the phase space [Fig. [3]Jc)] the alternate oscillations 
are a stable limit cycle C\ which surrounds the unsta- 
ble bidirectional state B\. When the bias current J is 
increased, two unidirectional solution U\ t 2 appear, and 
testability between the latter and C\ is possible (Region 
V) . This scenario corresponds to Fig. [3^d) . The basins of 
attractions of C/ lj2 are separated by the stable manifold 
of the saddle B\. A second unstable cycle C2 separates 
the basin of attraction of C\ from the basins of attraction 
of E/i 2. When noise is present in the system hopping be- 
tween U\ and U2 as well as hopping between L/1.2 and the 
alternate oscillations are possible, allowing the system to 
burst into periodic oscillations. A further increase of the 
current leads to the disappearance of C\ and the SRL op- 
erates in a bistable regime between the two unidirectional 
modes U\ t i (Region III). The phase space corresponding 
to this regime is shown in Fig.|3jb), similar to the case of 
4>k = 0.2. The stochastic terms induce hopping between 
Ui and U2, but no periodic oscillations can appear. 

The experiments have been performed on an InP- 
based multiquantum-well SRL with a racetrack geome- 
try and a free-spectral-range of 53.6 GHz. The exper- 
imental set-up is shown in Fig. [T] The device operates 
in a single-transverse, single- longitudinal mode regime at 
wavelength A = 1.56/im. The chip containing the SRL 
is mounted on a copper mount and thermally controlled 
by a Peltier element which is stabilized at a temperature 
of 28.57°C with an accuracy of 0.01°C. A bus waveg- 
uide made of the same active material as the ring has 
been integrated on the chip in order to couple power out 
from the ring. To this waveguide, an independent elec- 
trical contact has been applied. Sending current through 
the waveguide reduces the absorption. The power emit- 
ted from the chip is collected with a multimode fiber 
and detected with a 2.4GHz photodiode connected to an 
oscilloscope. The strength K and phase <f>k of the cou- 
pling between CW and CCW modes are not controllable 
during the fabrication process and they are a priori un- 
known. However, by using the cleaved facet of the fiber 
as a mirror, we are able to reflect power from one mode 
(for instance CCW) back into the waveguide and finally 
to the counterpropagating mode in the ring. The amount 
of power that is coupled to the CW mode can then be 
controlled by tuning the current I w on the waveguide, 
whereas its phase can be tuned by positioning the fiber 
facet with a piezoelectric controller. With this technique, 
we have achieved full control of the coupling parameter 
(j>k as well as the coupling strength K. 

We have fixed the waveguide current I w = 9.81mA and 
tuned the voltage on the piezoelectric controller until a 
symmetric operation of the SRL was achieved. Accord- 
ing to the discussion above, this corresponds to fixing 
the <pk parameter of the SRL. We now tune the bias cur- 
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FIG. 5. Experimentally measured time series for differ- 
ent bias currents I v and I w corresponding to the dynamical 
regimes of Fig. [3] (a) I w = 9.81mA, I v = 44.35mA; (b) I w 
= 9.81mA, I p = 45.21mA; (c) I w = 9.81mA, I p = 46.39mA; 
(d) I w = 12.0mA, I p = 40.23mA; (e) I w = 12.0mA, I p = 
44.34mA; (f) I w = 12.0mA, I p = 45.94mA. The insets show 
a zoom on relevant segments of the time series. 



rent I p on the ring, in order to reproduce the dynami- 
cal regimes predicted by the theory. We choose here to 
measure with the oscilloscope the emission in the CCW 
mode; due to the anti-correlated dynamics of the coun- 
terpropagating modes [T7], any change in the power in 
the CCW mode corresponds to an opposite change in 
the emission in the CW mode. A high (low) amplitude 
signal on the scope thus corresponds to operation in the 
CCW (CW) mode, whereas bidirectional operation ap- 
pears as an intermediate amplitude signal. Examples of 
timetraces of the CCW mode are shown in Fig. [5] (a)- (f). 
The ring reaches its lasing threshold at 34mA. For cur- 
rent values close to threshold, bidirectional lasing is ob- 
served (not shown). When increasing the current above 
a critical value I p ~ 44mA, hopping between the bidi- 
rectional regime and the two unidirectional modes ap- 
pear. Segments of time traces for I p = 44.35mA and 
Ip = 45.21mA are shown in Fig. [5|a)-(b). The average 
residence time in the CW and CCW state increases with 
the pump current, while the residence time in the bidirec- 
tional mode decreases. Testability between bidirectional 
and unidirectional modes is observed. In agreement with 
the phase space picture in Fig.|2ja), hopping events pref- 
erentially occur between the bidirectional and one of the 
unidirectional modes. A detail of a sequential transition 
from the CW mode to bidirectional operation, to the 
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CCW mode for I p = 45.21mA is shown in Fig. Wc). The 
agreement with the numerical simulations of Eqs. 
as shown in Fig. |4jb) is clear. The observation of some 
direct transitions between unidirectional modes suggests 
that the noise-induced diffusion length is not negligible 
when compared to the size of the basin of attraction of 
the bidirectional mode. When the pump current is in- 
creased to I p — 46.39mA as shown in Fig. |5jc), no bidi- 
rectional operation is observed, and direct transitions be- 
tween CW and CCW modes are possible as predicted by 
the phase space portrait in Fig.[3jb). Further increase in 
the bias current correspond to an increase in the average 
residence time in the two unidirectional modes. As such 
the whole bifurcation sequence as in Fig. [2ja) has been 
experimentally reconstructed. 

We then fix the waveguide current to I w = 12.0mA 
and we adjust the voltage on the piezoelectric until the 
symmetry in the system is restored. Once again we tune 
the pump current I p on the ring and we investigate the 
different dynamical regimes. Typical time traces for the 
CCW mode are shown in Fig. [5](cl)- (f ) . Close to thresh- 
old the laser operates in bidirectional regime, for slightly 
higher values of the bias current, alternate oscillations 
appear as reported in Fig. [5^d) for I p = 40.23mA. This 
operating regime reveals the phase space portrait shown 
in Fig.^c). When increasing the bias current, the am- 
plitude of the alternate oscillations increases until the 
SRL becomes multi-stable and the alternate oscillations 
coexist with the two unidirectional modes as shown in 
Fig. JBTc) for a bias current I p = 44.35mA. The inset in 
Fig. p[e) is a zoom on a burst that reveals the periodical 
oscillations. This operating regimes corresponds to the 
phase space portrait shown in Figj3^d). In this regime 
both transitions between the unidirectional modes and 
between the unidirectional modes and the limit cycle are 
observed, as allowed by the topology of the phase space 
in Figj^d). When the current is further increased, the 
bursts of oscillations disappear and bistability between 
CW and CCW modes is achieved as shown in Fig. ) 
for I p — 45.94mA. Such regimes corresponds to the phase 
space pictured in Figj3^b). Short excursions from CW 
to CCW operation and vice versa are observed in the 
time traces when the laser operates in this regime [see 
inset in Fig. [5jf )] . They have been previously observed 
[22j and explained as noise induced diffusion between the 
folds of the stable manifolds of the saddle point B\ in 
Fig. [3][b). The presence of such excursions in the time 
trace of Fig. [5jf ) represents a further confirmation of the 
phase space structure described in Fig. [3^b) . 

In conclusion, we have performed a comprehensive ex- 
perimental investigation of the phase space of a SRL. 
The dynamical regimes that we revealed experimentally 
match the phase space topologies of Eqs. (JTJ-([2|. A con- 
trol scheme based on the reflection of power from the 
cleaved facet of an optical fiber and the active bias of the 
bus waveguide has been devised in order to control the 



(otherwise unaccessible) coupling parameter </>/.. In this 
way, we could explore the dynamics of the SRL over the 
whole parameter space J—(f)k, including previously undis- 
closed regimes. Whereas alternate oscillations, bidirec- 
tional and unidirectional operation have previously been 
reported in SRLs [17], Fig.|5je) represents the first exper- 
imental observation of the coexistence between alternate 
oscillations and unidirectional operation, and more gen- 
erally shows coexistence of a limit cycle and two stable 
nodes in a semiconductor laser system. Moreover, coex- 
istence of three stable nodes has been demonstrated in 
Fig. [5^a)-(b).Our experiments have been performed on 
a SRL, and we expect our findings to be extendable to 
other systems that display Z2 symmetry, such as for in- 
stance disk-lasers. 
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